Electronic properties of guanine-based 



nanowires 



A. Calzolari, * R. Di Felice, and E. Molinari 

INFM-S 3 - National Research Center on nanoStructures and Biosystems at 
Surfaces, and Dipartimento di Fisica Universitd di Modena e Reggio Emilia, 

HI 100 Modena, Italy 



Abstract 

We present a first-principle study of the electronic and conduction properties of 
a few classes of nanowires constituted of guanine (G) molecules, self- assembled in 
different geometries. We first analyze the effect of the vertical tt-tt interaction in 
model G-stack columns. Then, we exploit the results obtained from those models to 
interpret the features of realistic stacked and hydrogen-bonded structures, namely 
the guanine quadruple helices and the planar ribbons. With respect to natural DNA, 
the different structures as well as the inclusion of metal cations, drastically affect 
the bonding pattern among the bases, introducing novel features in the electronic 
properties of the systems. These supramolecular G-aggregates, alternative to DNA, 
are expected to show intersting properties for molecular electronics applications. 
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1 Introduction 



The idea that single molecules and/or molecular arrays might be anchored to 
inorganic supports to be integrated into electronic circuits has inspired many 
scientists in the last decade. An intense research effort is being developed to 
identify the suitable molecular candidates, exploit them to fabricate devices, 
and test the device efficiency. This is the basic task of molecular electronics, 
which is the solution of choice proposed for the continued miniaturization 
of electronic circuits beyond the downfall of the Moore's law predicted for 
silicon-based electronics. 

By virtue of their recognition and self-assembling properties, [1] DNA molecules 
appear as particularly suitable candidates as components for electronic devices 
at the nanoscale. DNA's natural function is information storage and transmis- 
sion, through the pairing and stacking characteristics of its constituent bases. 
The idea that it may also carry an electrical signal would be extremely inter- 
esting in the bio-chemical environment, but would also be appealing for the 
exploitation of DNA in the nanotechnologies. 

In 1962, D.D. Eley and D.I. Spivey [2] suggested that tt — tt interactions be- 
tween stacked base pairs in double-stranded DNA could provide a pathway 
for rapid, one-dimensional charge separation. The issue of charge migration in 
DNA has recently become a hot topic in solution chemistry experiments. [3] 
Most of the available experimental results concern the migration of positive 
charge carriers (e.g. radical cations) within DNA molecules in solution. A 
major target for oxidation is guanine (G), [4,5] the natural nucleic-acid base 
(Fig. la) with the lowest ionization potential. [6] The discussion of the results 
obtained using different methodologies and systems focused on the depen- 
dence of the charge-transfer efficiency on the length of a 7r-pathway, serving 
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as a bridge between remote donor and final acceptor sites. In very general 
terms, two different regimes can operate in charge conduction through DNA 
molecules: as a function of the base sequence, [7,8] both short-range (a few 
A) and long-range (up to ~ 200 A) charge migration have been observed. In 
particular, in guanine-rich sequences of DNA, the oxidation/trapping process 
can happen over a long distance between the donor and the acceptor sites. 

On the other hand, the construction of solid-state DNA-hybrid devices is a 
very recent achievement. Both the combinatory principles intrinsic to nucleic 
acids and their chemistry can be exploited in building precise, miniaturized 
and locally modulated patterns, where the drawing of functional arrays is 
obtained through a series of programmed chemical reactions rather than by 
physical handling of samples. The realization of wired architectures via self- 
assembly opens the way to scale the size of the devices down to the molecular 
level. [9] In principle, the DNA-based materials may be used either as conduc- 
tive wires or as templates for other materials. [10] While important successes 
have been achieved in the last years in the immobilization of selected DNA 
molecules onto substrates and between nanoelectrodes [11] confirming their 
excellent mechanical properties, a clear understanding of DNA conductivity 
and its microscopic origin has not yet been achieved. A recent series of ex- 
periments, which provided the direct measurements of the dc conductivity in 
DNA-hybrid structure, reports contrasting results: [12,13] depending on base 
sequence, molecule length, environmental conditions, substrates, and electrode 
materials, DNA exhibits insulating character, [10,14,15] semiconductor-like 
transport characteristics, [16,17] ohmic behavior, [18] and proximity-induced 
superconductivity. [19] 

Given these contradicting reports, no definite conclusions on the intrinsic con- 
duction properties of DNA can be presently drawn from solid state experi- 
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ments. As found for the charge transfer rates measured in solution chemistry 
experiments, the electrical properties of DNA-based devices are strongly de- 
pendent on the sequence of the base pairs: in particular an abundance of gua- 
nines seems to enhance the conductivity of the double helix [16,20]. However, 
even in the cases in which charge transport is observed, the current flow is very 
low, so that DNA might not to be reasonably considered a good conducting 
wire. 

In this paper, we present a few nucleobase-aggregates, alternative to natural 
DNA, expected to show interesting properties for molecular electronics ap- 
plications. Based on the experimental results stated above, we focus on well- 
defined molecular structures by restricting the base sequence to the guanine 
base alone. Self-recognition and self-assembly processes of G-aggregates are 
controlled by a unique sequence of functional groups, which act as donors and 
acceptors of H-bonds (Fig. la). Among the bases of DNA, guanine undergoes 
extensive self-assembly, mediated by H-bonding between G's, to give different 
aggregation states (dimers, quadruple helices, ribbons) depending on the sol- 
vent concentration and on the presence (or absence) of specific cations. [21,22] 
With respect to native DNA, the different structures drastically affect the 
bonding pattern among the bases, introducing novel features in the electronic 
properties of the system. 

In order to understand the complicated electronic properties of these mate- 
rials, we first focus on the effects of the n — n coupling in selected model 
systems, constituted of vertical stack of guanine (poly(G) columns). Then, 
we consider two realistic systems: the stacked G-quartet helices (G4~wires or 
quadruplexes) [22,23] and the planar H-bonding G-ribbons. [24] The latter are 
the building blocks of guanosine fibers, successfully exploited in the construc- 
tion of nanodevices. [25,26,27] 



4 



By means of ab-initio calculations, we study the electronic and conduction 
properties of these biomaterials from a microscopic point of view. We focus 
on the effects of the it — ix coupling and of the hydrogen bonding on the 
electrical properties of the self-assembled guanine nanowires. In particular, 
the question that we address is whether the electronic properties of extended 
G-based structures can be described in terms of extended orbitals and account 
for a band-like mechanism for charge transport. 

2 Method 

In our first principle study of the self-assembled guanine nanowires, the elec- 
tronic structure is described according to the Density Functional Theory (DFT) [28]. 
DFT is lately gaining a large credit in the scientific community as a reliable 
and accurate method to describe large-scale biomolecular aggregates [29] in- 
cluding guanine-based stacks [30,31,32] and native DNA [14,33,34]. We include 
the Generalized Gradient Approximations (GGAs) to the exchange-correlation 
(XC) functional, in order to improve the description of highly directional in- 
teractions, such as H-bonding in molecular systems. [35] 

The electron-ion interactions are simulated via ab-initio pseudopotentials for 
each chemical species (H, C, N, O, K). The electron wavefunctions are ex- 
panded in a basis of plane waves. We present two set of calculations with 
slightly different technical details: poly(G) columns and ribbons are treated us- 
ing BLYP [36] XC-functional, norm-conserving (Troullier-Martins [37]) pseu- 
dopotentials in the in the Kleinman and Bylander factorized form [38], and 
an energy cut-off of 50 Ry for the expansion of the single-particle wavefunc- 
tions. In G4-wire calculations, instead, we exploit PW91 [39] XC-functional, 
ultrasoft [40] pseudopotentials, and 25 Ry energy cut-off for the wavefunction 
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expansion. Initial test on gas-phase molecules and simple G-based aggregates 
confirmed that the latter approach does not modify substantially the descrip- 
tion of the electronic structure of the systems, while allowing the treatment 
of larger unit cells. 

We simulate one-dimensional (ID) crystals by periodically replicating guanine 
building-block aggregates along the axis of the nanowire, and including a thick 
layer of vacuum in the other spatial directions to avoid spurious interactions 
between adjacent replicas of the system. 

For each selected G-assembly, we optimize the atomic geometry and the elec- 
tronic structure. We perform quenched ab-initio molecular dynamics simula- 
tions, according to a total-energy-and-force scheme for the search of the local 
minima. All the atoms are allowed to relax, until the forces vanish within an 
accuracy of 0.05 eV/A. For each structure, we obtain both the geometry and 
the corresponding single-particle electron energies and wavefunctions. 

Based on previous experimental and theoretical results [3,41] about similar 
nucleoside aggregates, indicating that the sugar-phophate backbone is not di- 
rectly involved in the electron transport processes, we focus on pure G aggre- 
gates that, with respect to the experimental structures, neglect the external 
backbone (G4-wire), the sugar and the alkyl chains (G-ribbons), and the sur- 
rounding medium. The substituent groups attached to guanines are replaced 
with H atoms (H# in Fig. la) to saturate the dangling bonds. 

3 Poly(G) model columns 

In order to clarify the complex intramolecular interactions in nucleobase- 
aggregates, we started by investigating the role of n-n coupling in simple 
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toy systems. [30] We considered several infinite periodic stacks of guanines in 
a wire-like assembly, with different inter-plane configurations. We stress that 
these poly(G) columns do not describe real systems existing in nature, but they 
are useful tools to understand the effects of 7r interactions as a function of the 
geometrical parameters. The results obtained with these theoretical models 
constitute a first step towards the interpretation of the electronic features of 
some real systems, e.g. guanosine fibers and G-quartet nanowires. 

We started our simulation on isolated stacked dimers, pairs of G's lying in 
parallel planes whose distance in the perpendicular direction is an output 
of our calculations. We analyzed several configurations, characterized by the 
relative azimuthal rotation angle of the two G's in the pair, around the axis 
identified by the stacking direction, and/or by an in-plane translation. The 
different geometries realize different tt—tt overlaps. After relaxing the structure 
of the G-dimers, we determined the equilibrium interplanar distance to be 
used for periodic column calculations. For all the computed configurations 
the vertical distance between consecutive stacked planes is 3.37 A, which is 
the experimental distance observed in the B-DNA and in the G-quadruplexes. 
The detailed description of all selected configurations and their properties are 
reported elsewhere. [30] In the following, we present the results for two extreme 
cases with a very large and a very small % — ir superposition. 

In the first configuration (Fig. 2, left inset) the two G's of the repeated dimer 
are perfectly eclipsed; this realizes the maximum of the n — n superposition. In 
the second one (Fig. 2, right inset), the relative rotation angle is 36° around a 
suitable axis so to mimick the typical configuration of two consecutive guanines 
along a double strand of B-DNA. In this case, the n — n overlap is very small. 

The analysis of the electron properties reveals interesting features. We cal- 
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culated the band structure along the stacking direction for both geometries 
reported in Figure 2. [42] In both cases, we obtained a wide energy band-gap 
- strongly underestimated in DFT - between the valence band maximum and 
the conduction band minimum. 

In the eclipsed configuration (maximum n — n overlap), the bands deriving 
from the Highest Occupied Molecular Orbital (HOMO) and from the Lowest 
Unoccupied Molecular Orbital (LUMO) are dispersive. The calculated effec- 
tive masses for holes and electrons (m^/m =1.04, m e /m =1.41) are similar to 
those of inorganic wide band-gap semiconductors, such as group-Ill nitrides, 
for which a band-like conductivity is demonstrated (in presence of doping). 
The n — n interactions along the stacking direction induce extended channels 
for charge migration through a band transport mechanism. This interpreta- 
tion is in agreement with the character of electronic states near the energy 
gap: In the columnar arrangement, through the interplanar G-G coupling, the 
7r-like HOMO of guanine (Figure lb) induces the formation of Bloch-type de- 
localized orbitals that extend along the stacking directions (Figure 3). Thus, 
we conclude that the electronic properties of the eclipsed G-wire may sup- 
port band-like electron transport through the base stack (provided a suitable 
doping mechanism might be achieved). 

The picture drastically changes if we consider the B-DNA type configuration 
(small 7r — 7r overlap) . The energy bands are dispersionless (right panel of Fig. 
2) and the effective masses are very high (m/ l /m =5.25, m e /m — > oo). The 
low 7T — 7T superposition does not induce the formation of extended Bloch-type 
orbitals along the stacking direction, and the electronic states are localized 
at the individual molecules. This configuration does not support a band-like 
conductivity, and has electronic properties typical of a molecular insulator. 
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In synthesis, the electronic properties of vertical stacks of G molecules are 
very sensitive to n — n coupling: a large t\ — tx superposition is a necessary 
condition to sustain a Bloch-like conduction along the stacking direction. 

Although these poly(G) structures do not directly simulate any real molecule, 
they provide some general notions that are helpful to interpret experiments on 
G-rich DNA molecules. For instance, our results suggest that relative geom- 
etry of neighboring bases in B-DNA does not support energy dispersion and 
orbital derealization (Bloch-type conduction) in a native poly(dG)-poly(dC) 
oligomer, [16] due to the small ti-ti overlap. However, one cannot exclude that 
the interactions with the backbone and the external environment, as well as 
the thermal fluctuations, may locally change the base-pair superposition for a 
finite time, thus favoring an increase of the n — n overlap and the formation 
of extended orbitals delocalized over a few bases, and eventually generating 
polarons. 

As a final note to the above discussion, we can state that band-like conductiv- 
ity (due to dispersive energy bands) can only give a small contribution [41] (if 
any) to the phenomenon of long-range charge migration in DNA. Other rele- 
vant concurring mechanisms (e.g. tunneling, thermal hopping, polaron drift) 
must be taken into account to obtain a comprehensive interpretation. 

4 G4-wires 

The poor conductive properties of DNA (as suggested by the experiments [12] 
and by the results of previous section) have driven the research towards novel 
biomaterials that might provide electrical properties suitable for molecular 
electronics, while maintaining the addressability and structured arrangement 
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that are so attractive in the native DNA. We identify a possible candidate 
in a homoguanine DNA-derivative: the G4-wire (or quadruplex). Due to the 
low ionization potential of guanine, this form might be suitable to mediate 
charge transport by hole conduction along the helix. This low-dimensional 
structure, widely studied for its biological properties, [43] exhibits a surprising 
high stability in different chemical conditions, and may form nanoscale (100- 
1000 nm long) ordered wires. [44] 

In the presence of appropriate metal cations (especially K + and Na + ) solutions 
of homoguanylic strands, and of some lipophilic guanosine derivatives self- 
assemble in right-handed quadruple helices (Fig. 4a), constituted of stacked 
supramolecular tetramers, known as G-quartets (G4). Each G4 (Fig. 4b) con- 
sists of four coplanar hydrogen-bonded guanines, arranged in a square-like 
configuration. The metal cations, hosted in the central cavity resulting from 
the stacking of the G-quartets (Fig. 4c), play a leading role in the assembly and 
stability of the aggregate. They also call for the investigation of the effects of 
the metal-molecule interactions on the electrical properties of the nanowires: 
such interactions are currently explored as a viable route to enhance the con- 
ductivity of DNA derivatives. [17] 

The starting atomic configuration for our simulations was extracted from the 
X-ray analysis of the short quadruplex d(TG<iT). [23] In the definition of 
the unitary cell, we exploited the periodicity modulo three of the G-quartet 
stack (Fig. 4c): we simulated periodically repeated G4- wires by means of three 
stacked tetrads and three intercalated K + cations as shown in Figures 4a, c. 
Details about the geometry and the thermodinamic stability of G4-wires are 
reported elsewhere. [31] 

We analyze of the electronic structure in terms of the molecule-molecule and 
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the metal-molecule interaction. In the description of the guanine-guanine cou- 
pling we can also distinguish between the in-plane and the out-of-plane inter- 
actions. Indeed, the in-plane structure of G-quartet (Fig. 4b) introduces novel 
features related to presence of the H-bonds [30,31], while the out-of-plane be- 
haviour is controlled by the tt-tt coupling along the stacking direction. The 
results of the simulation show that H-bonding is not capable of modifying the 
internal configuration of the single guanine; the polar nature of the H-bonds 
does not involve the sharing of electronic charge as the covalent bonds do. The 
molecular orbitals of the guanines do not laterally couple to form an in-plane 
supramolecular orbitals for the whole tetramer; the same behavior is observed 
in the isolated G4, as illustrated by the HOMO state plotted in Fig. 4d. We 
conclude that the planar H-bonds do not significantly affect the band-like 
conductivity of the whole helix. 

In the previous section, we showed that tt-tt coupling may give rise to de- 
localized Bloch-type orbitals, whose band dispersion depends on the relative 
rotation angle between nucleobases in adjacent planes. In the guanine quadru- 
plexes, the square symmetry of the tetrads increases the spatial tt-tt overlap 
with respect to a segment of G-rich B-DNA (30° instead of 36°), and one 
might hope to find an enhanced band-like behavior. Instead, we find that the 
out-of-plane tt superposition is not sufficient to induce the formation of de- 
localized orbitals along the axis of the quadruplex. Consequently, the energy 
bands along the stacking direction are flat, [32] reflecting the fact that the 
orbitals remain localized at the component molecules. 

However, the band structure displays the presence of manifolds, due to the 
intermolecular interactions that split the energy levels. The electron orbitals 
associated to a manifold have identical character and localize on the con- 
stituent G molecules of the cell. The energy levels in a multiplet are separated 
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by an average energy difference of 20 meV; therefore, even the coupling with 
a weak external perturbation (e.g. thermal fluctuations, electric field, etc.) 
might be sufficient to allow the G-localized orbitals of a multiplet to inter- 
act, producing an effective delocalized orbital. This interpretation suggests a 
description of the electronic structure where the spreading of the energy lev- 
els of the manifolds leads to the formation of effective band-like peaks in the 
Density of States (DOS) (Fig. 5, black area). The DOS of the G4-wire thus 
appears similar to those of wide-bandgap materials with dispersive band-like 
structures. We conclude that the coupling of the electronic states belonging 
to a manifold is coherent with the formation of effective extended electron 
channels, suitable to host mobile charge carriers along the wire. [32,45] 

The comparison with the empty G4 tube (without K + ions) [31] and with 
DNA structures available in the literature [14,35] underlines how the forma- 
tion of dense manifolds (band-like peaks) is a general feature of the base-base 
interactions in the H-bonded nucleobase stacks. Instead, the details of inter- 
planar coupling depend on the specific arrangement (base sequence, rotation 
angle) and are influenced by the presence of metal cations intercalating the 
base stack. [32] Additionally, the effects due to the presence of the cations 
depend on the specific choice of the metal species. 

In the case of potassium, the detailed description of its electronic structure 
(which includes the explicit treatment of 3p semicore shell) allows to in- 
vestigate the metal-molecule interactions, neglected so far. We analyzed the 
metal-base coupling, projecting the total density of states of the quadruplex 
(black area in Fig. 5) on guanine (red) and potassium (green) atomic orbitals. 
Whereas for a large part of the spectrum the whole contribution to the DOS 
stems only from the guanine states, the inclusion of K + ions strongly affects 
the electronic properties in the energy region near the Fermi level (E=0 in 
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Fig. 5), as indicated by the superposition of G and K contributions in Figure 

5. Potassium contributes to the G4-wire electronic structure with partially 
occupied sp-orbitals, which derive from the mixing of the atomic 4s and 3p 
states. The ionization process uniformly decreases of one unit the total num- 
ber of electrons, without changing the character of these mixed states, which 
further couple with the HOMO's of constituent guanines giving rise to the 
effective-HOMO of the entire helix (marked with an arrow in Fig. 5). The 
interaction between hybridized sp-orbitals of K + and the guanine orbitals ac- 
complishes a double effect. First, with respect to the empty G4 aggregate, 
the inclusion of potassium enhances the conduction properties of the system 
increasing the HOMO band width, and generating further extended electron 
channels available for charge (hole) transport through the nanowire. Second, 
the rearrangement of the mixed K-G orbitals gives origin to a partially filled 
HOMO-derived peak, that characterizes the G4-wire as an intrinsic p-doped 
material. The observation of delocalized holes at the topmost valence band 
is in agreement with the description of the electronic structure of potassium 
discussed above, where each K + ion contributes to the whole electronic struc- 
ture not with three fully occupied p-levels, but with four partially occupied 
sp-states. 

We summarize the complex electronic properties the quadruplexes in Figure 

6, where we report the contour plot of the HOMO-derived state of the wire 
in two different planes. The left panel of Fig. 6 has been obtained by slicing 
an isosurface of the HOMO charge with a plane perpendicular to the axis of 
the helix and containing a G-quartet. The maxima of charge density (orange 
color) are localized around each guanine and no supramolecular orbitals are 
visible across the H-bonds. The in-plane charge distribution is basically the 
same of the isolated G-quartet plotted in Figure 4d. The right panel of Figure 
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6 shows, instead, the contour plot of the HOMO-derived state in a plane 
parallel to the helix and containing the central metal cations. We recognize a 
uniform distribution of the n orbitals deriving from the single G's, but also the 
inner states stemming from the K + ions. Two types of delocalized (Bloch-like) 
spatial electron channels can be identified. The first one, which is extended 
through the guanine core of the helix, is due to the molecule-molecule stacking 
interaction and corresponds to the formation of the effective band-peaks in the 
density of states. Similar channels are also observed for the other manifold- 
derived peaks. [32] Conversely, the second type of channel is specific of the 
HOMO peak, because it results from the metal-molecule interaction, which 
is active only in the energy region near the top of valence band. This hybrid 
orbital is centered around the potassium ions in the central cavity of the wire, 
and it clearly shows the K + coupling with the coordinated oxygen atoms of 
the G molecules. 

The above results, along with their attractive mechanical and self-assembling 
properties, suggest that the G4- wires should be explored as viable DNA-based 
conductors for nanoscale molecular electronics. 

5 G-ribbons 

As mentioned above, guanine derivatives may follow different self-assembly 
schemes, mediated by H-bonding. Besides the formation of the Hoogsteen 
dimers and four stranded helices, a planar phase has been observed: in con- 
trolled experimental conditions, lipophilic guanosine derivatives self-assemble 
in a ribbon-like structure, both in solution and in the solid-state fibers [21]. 
The guanosine fibers [46] have been recently used to fabricate working nan- 
odevices such as rectyfing diodes [26] and transistors. [27] 



14 



The planar ribbons are self-assembled quasi-lD structures of G's linked by a 
H-bonding network. We firstly studied the isolated planar ribbon, [48] whose 
structure is presented in Figure 7. To simulate extended ribbons we assem- 
bled coplanar G molecules according to experimental geometry [24]. In the 
same way as we did for columnar G-stacks, we obtained the intermolecular 
distance by optimizing the structure of an isolated dimer. We used this value 
to fix the periodicity of the extended structures, whose geometry was also op- 
timized. The periodicity of the ribbons (l.lnm) is in good agreement with the 
experimental data (1.2nm). [24] The donor-acceptor distances (NH- • • N=2.9 
A, NH- • • 0=2.8 A) are consistent with previous theoretical results for sim- 
ilar systems (e.g. DNA base-pairs) [47]. We find that the main effect of the 
H-bonding is to stabilize the structure, as proved by a large calculated energy 
gain (A£ /orm =-820 meV/G). 

The analysis of the electronic structure shows that, as it happens for the G- 
quartet, the H-bonds do not favor the formation of supramolecular orbitals 
extended on the whole ribbon. On the contrary, electronic states remain lo- 
calized around single G's. Such electronic structure is not compatible with a 
Bloch-like scheme for electron transport: in fact the energy bands along the 
axis of the ribbon are totally flat, and the calculated density of state (not 
shown) [48] is characterized by sharp peaks corresponding to energy levels 
of the single guanines. We conclude that an isolated planar ribbon should 
behave as an insulator, [30,48] unless charge carriers migrate according to 
discrete hopping process. 

A mechanism concurring to charge transport may derive from the dipole-dipole 
interaction, which arises along the ribbon. [48] Indeed, G-ribbons have a net 
dipole oriented along their axes (Fig. 7); it derives from the arrangement of the 
dipole moments of each guanine, [49] which vectorially add up to yield a finite 
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and large spontaneous dipole moment lying parallel to the axis of the ribbon. 
The associated electric field may favour the charge motion through adjacent 
G's, by polarizing the medium. It is worth noting that the presence of the 
macroscopic dipole is a consequence of the particular guanine arrangement 
in the ribbon geometry. For instance, the square-like symmetry of the planar 
G-quartet (Fig. 4b) does not lead to a finite dipole for the whole tetrad. 

Finally, to gain further insight into the fundamental interactions in extended 
gaunine aggregates, we studied a 2D crystal of stacked ribbons. Following the 
results of Sec. 3, we expect that the electronic properties of such systems vary 
as a function of the relative base-base superposition. Here, we present the 
results for a structural model in which the ribbons are perfectly eclipsed and 
separated by 3.4Aalong the stacking direction. In principle, this geometry does 
not reproduce the experimental solid-state arrangement of the ribbons: Even 
though experimental results prove that the ribbons self-organize in fibers, the 
details of the stacking arrangement may depend on the length of the attached 
alkylic chains. Thus, while our eclipsed model does not allow for a quantitative 
evaluation of the conduction properties of the real system, it allows us to 
underline the relative role of the H-bonding and the % — n interactions when 
the n — tt coupling is maximized. 

The electronic properties of the stacked configurations are simply the combi- 
nation of two interactions: the H-bonding on the ribbon plane and the % — n 
along the stacking direction. The band structure (middle panel of Fig. 8) 
shows two different behaviors along the ribbon axis (T — X in the picture) and 
along the stacking direction (F — A). In particular, as for the isolated ribbon, 
energy bands are almost flat along the axis (T — X); instead, energy disper- 
sion occurs along the stacking direction (T — A), as it happens for eclipsed 
poly(G) columns. This is consistent with the analysis of single-particle orbitals 
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shown in Figure 8. Whereas H-bonding does not favor lateral coupling, and 
the electronic states remain localized at the G's (top panel of Fig. 8), the 
n — n coupling induces the formation of Bloch-type orbitals delocalized along 
the vertical direction (bottom panel of Fig. 8). 

Stacked eclipsed ribbons present a wide-bandgap semiconducting behavior, 
compatible with a band-like conduction in the stacking direction. As it hap- 
pens with vertical columns, the band properties of stacked ribbons depend on 
the amount of the ty — tv overlap: different geometries with a lower n — n super- 
position, [30] are expected to have worse band-like properties with larger effec- 
tive masses with respect to the eclipsed configuration. However, as a function 
of the superposition, the ty — it interaction may enhance the band properties, 
improving the mobility of charge carriers only along the stacking direction. 

In summary, H-bonding and 7r— 7t coupling play different roles: H-bonding con- 
trols the in-plane self-assembly of guanines and stabilizes the ordered struc- 
ture, whereas n — tt coupling has a stronger influence on the out-of-plane 
electronic properties and determines the contribution of delocalized band-like 
orbitals to the overall charge transport of G-rich systems. 

Our results account for the experimental evidence of a semiconducting behav- 
ior of ribbon-like fibers, [26] where even a partial guanine overlap may favor 
the formation of band-like channels available for charge transport. 

6 Conclusions 

We presented a microscopic theoretical investigation of the electronic proper- 
ties of a few selected nanowires obtained by the different spatial arrangements 
of guanine molecules. 
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By an extensive study of guanine model solids, we examined the characteristics 
of 7r-7r coupling and H-bonding, and discussed their effect on the conductiv- 
ity of the molecular assemblies. We found that Hydrogen-bonding does not 
support the band transport: no band dispersion is present along H-bonds in 
planar guanine aggregates, such as planar ribbons or isolated G-quartets. In- 
stead, base stacking is accompanied by n-n interactions that, for sufficiently 
large overlap between adjacent tt orbitals, may induce energy dispersion and 
are consistent with a finite contribution to coherent transport in the pres- 
ence of doping or photoexcitation. Therefore, band transport may be par- 
tially responsible for charge mobility in nucleotide aggregates, in structures 
characterized by a large base-base superposition: the electronic properties of 
such materials are similar to those of inorganic wide-bandgap semiconductors. 
These basic features occur in all the different homo-guanine structures that 
we have considered (poly(G) columns, G-ribbons, G4- wires). However, the de- 
tails of the electronic properties in each nanowire are strongly affected by the 
specific supramolecular aggregation state and the coupling mechanisms such 
as the metal-molecule hybridization in the quadruple helices or the dipolar 
interaction in G-ribbons. 

These features, along with the possibility of forming extended stacked wires 
at the nanoscale length, make both the G4-wires and the G-ribbons appeal- 
ing materials for the development of biomolecular electronics, possibly more 
promising than DNA. 
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Fig. 1. Isolated guanine molecule, (a) Atomic structure; (b) isosurface plot of HOMO 
state. Different colors stand for different chemical species. The hydrogen Hr replaces 
the substituent groups attached to the molecule in the experimantal systems. 
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Fig. 2. Band structure for eclipsed (left) and B-DNA type (right) geometries, cal- 
culated along the stacking direction. The single particle energies reported in these 
plots are relative to the top of the highest valence band. The insets show the top 
view of the two configurations respectively. 
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Fig. 4. (a) Side view of periodically repeated G4-wire; (b) Geometry of isolated 
G-quartet; (c) Schematic description of the periodicity modulo three in G-quartet 
stacks; (d) Isosurface of the HOMO state of the isolated G-quartet. Colors refers to 
different chemical species as in Fig. 1, green stands for K ions. 
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Energy (eV) 

Fig. 5. Total and Projected Density of States of the wire. The shaded black area 
represents the total DOS. The squared red (dotted green) area marks the projection 
of the DOS onto G (K) contribution. Note that the total DOS and the G-DOS 
coincide in a large part of the spectrum, where the red area almost coincides with 
the black one. The zero of the reference energy scale indicates the Fermi level. 
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Fig. 6. Contour plot of the symmetrized charge density of the states belonging to the 
HOMO manifold, in a plane perpendicular to the axis of the helix and containing a 
G-quartet (left); and parallel to the helix and containing the central metal cations 
(right). 
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Fig. 7. Isolated planar G-ribbon: relaxed geometry, and dipole moment (p). 
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Fig. 8. Central panel): Band structure of eclipsed stacked ribbons along the axis 
(r — X) and the stacking directions (r — ^4). The zero energy is fixed to the top of 
valence band. Isosurface of HOMO state in the plane containing the ribbon (top), 
and in a plane perpendicular to the ribbons and parallel to the stacking direction 
(bottom). 
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